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We show how to build a multi-scale entanglement renormalization ansatz (MERA) representa¬ 
tion of the ground state of a many-body Hamiltonian H by applying the recently proposed tensor 
network renormalization (TNR) [G. Evenbly and G. Vidal, arXiv: 1412.0732 to the Euclidean time 
evolution operator e~^ H for infinite (3. This approach bypasses the costly energy minimization of 
previous MERA algorithms and, when applied to finite inverse temperature fi, produces a MERA 
representation of a thermal Gibbs state. Our construction endows TNR with a renormalization 
group flow in the space of wave-functions and Hamiltonians (and not merely in the more abstract 
space of tensors) and extends the MERA formalism to classical statistical systems. 
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Consider a strongly interacting quantum many-body 
system in D spatial dimensions described by a micro¬ 
scopic Hamiltonian H. Understanding its collective, low 
energy behavior is a main goal in condensed matter and 
high energy physics, one that poses a formidable theoret¬ 
ical challenge. To tackle this problem, plenty of methods 
have been proposed based on the renormalization group 
(RG) EH , that is, on studying how the physics depends 
on the scale of observation. Weakly interacting systems 
can be addressed perturbatively using momentum space 
RG [2]. Instead, strongly interacting systems often re¬ 
quire non-perturbative, real space RG methods, as pio¬ 
neered by Kadanoff and Wilson Q. 

Improving on Kadanoff and Wilson’s proposals, 
White’s density matrix renormalization group (DMRG) 
for quantum spin chains Q established how to systemati¬ 
cally preserve the ground state wave-function during real- 
space coarse-graining - namely, by preserving the sup¬ 
port of its reduced density matrix. Similarly, Levin and 
Nave’s tensor renormalization group (TRG) [HJ taught 
us how to coarse-grain Euclidean path integrals of one¬ 
dimensional quantum systems (also partition functions 
of two-dimensional classical systems). Both DMRG and 
TRG are very successful, versatile approaches. However, 
they depart significantly from the spirit of the RG, in 
that they produce a coarse-grained, effective description 
of the system that still retains some irrelevant micro¬ 
scopic details. As a result, these methods (i) fail to define 
a proper RG flow, one with e.g. the correct structure of 
fixed points; and (ii) struggle to deal with critical systems 
or systems in D > 2 dimensions, where the accumulation 
of irrelevant microscopic degrees of freedom is more sig¬ 
nificant and harmful. 

These difficulties have been solved with two closely re¬ 
lated proposals. First, entanglement renormalization was 
put forward to address the above two problems in the 
context of ground state wave-functions [6, 7]. By intro¬ 
ducing disentanglers, which remove short-range entangle¬ 


ment, a proper RG flow is generated, as well as an RG 
transformation that is computationally sustainable even 
at criticality. In addition, entanglement renormalization 
leads to an efficient tensor network description of ground 
states for critical systems, the multi-scale entanglement 
renormalization ansatz (MERA) Ml, of interest not only 
as a many-body variational state [7- 1J] but also as a lat¬ 
tice realization of the holographic principle 

Tensor network renormalization (TNR), on the other 
hand, was more recently proposed to tackle the same 
problems in the context of Euclidean path integrals (and 
classical statistical systems) 0. The Euclidean path 
integral Z & tr ( e~^ H ^ is represented by a tensor net¬ 
work, consisting of copies of a single tensor A [l8[ , which 
extends both in space and Euclidean time directions. 
Through local manipulation of this tensor network, TNR 
produces a sequence of tensors, 

A A' A" -> A fp , (1) 


corresponding to increasing length scales, which flow to¬ 
wards some infrared fixed-point tensor Af p . The later 
retains only the universal features of the phase or phase 
transition [17]. Once again, the key of the approach is 
the removal of short-range correlations by disentanglers. 

In this Letter we establish a close connection between 
the two approaches. We show that, when applied to the 
Euclidean path integral restricted to the upper half plane, 
TNR generates a MERA for the ground state of Hamilto¬ 
nian H. More generally, TNR also produces a MERA for 
the thermal Gibbs state pg = e~^ H /Z at finite inverse 
temperature /3, as well as for the low energy eigenstates 
of H on a finite periodic chain. Our result provides an 
alternative route to the MERA, one that bypasses the 
costly energy minimization of previous algorithms fsl, 9j 
and has several other significant advantages, both con¬ 
ceptual and computational, that we also discuss. Among 
those, we emphasize that (i) we obtain the first correct 
MERA representation of the thermal state pg 0 , to- 
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gether with an algorithm to find it; (ii) the connection 
implies that TNR produces an RG flow in the space of 
wave-functions and Hamiltonians, and not just in the ab¬ 
stract space of tensors, Eq. [TJ (iii) the MERA formalism 
can be applied now also to classical statistical systems. 
For simplicity, we consider a translation invariant system 
in D = 1 dimensions, although the key results generalize 
to inhomogeneous systems in dimension D > 1. 

Tensor network for Euclidean time evolution. — Given 
a translation invariant Hamiltonian H in one dimen¬ 
sion, we use a standard procedure [Jo] to produce a 
two-dimensional tensor network representation of the Eu¬ 
clidean time operator e~^ H or Euclidean path integral 
tr (e~P H ). This tensor network is made of copies of a 
single tensor A. If both the system size L and the in¬ 
verse temperature (3 are infinite, then the network spans 
the entire (x, r)-plane, where x and r label space and Eu¬ 
clidean time, respectively. Here we will consider tensor 
networks for e~^ H on three different geometries, obtained 
by introducing a horizontal cut at r = 0 and by choosing 
L and [3 to be either finite or infinite 0. 

TNR yields MERA. — Let us start with the upper half 
plane, Fig. HJ which corresponds to the ground state |T) 
of Et on an infinite lattice. The network has an open 
boundary at r = 0, with an infinite row of open indices, 
one for each site of the one-dimensional lattice on which 
H acts. We first apply TNR everywhere on the upper half 
plane except near r = 0, where we keep the open indices 
of the tensor network untouched. TNR acts through an 
intrincate sequence of local replacements 0 . Here we 
skip the technical details [22] and focus instead on de¬ 
scribing the final result: a coarse-grained tensor network 
with effective tensor A! for most of the upper half plane, 
in accordance with Eq. HJ together with a double row 
of special tensors, so-called disentanglers and isometries, 
which correspond to one layer of the MERA [23]. These 
tensors connect the microscopic degrees of freedom at 
length scale s = 0 (represented by the original open in¬ 
dices of the network) with the coarse-grained degrees of 
freedom at length scale 5 = 1 (represented by the lower 
indices of the lowest row of tensors A'). We can now re¬ 
peat the process on tensors A', to obtain coarse-grained 
tensors A" and a second row of disentanglers and isome¬ 
tries, i.e. a second layer of the MERA, connecting scales 
5 = 1 and 2. Iteration then produces a full MERA ap¬ 
proximation for the ground state |T) of iL, encompassing 
all length scales 5 = 0,1, 2, • • •. 

Thermal MERA. — Let us now consider a horizontal 
strip of finite width /?, Fig. [2fa), which is proportional 
to the thermal state pp = e~^ H /Z. This time we have 
two boundaries, each with an infinite row of open in¬ 
dices: the incoming and outgoing indices of the Euclidean 
time evolution operator e~^ H . As before, we use TNR 
to coarse-grain the tensor network, except near its open 
boundaries, where we do not touch the open indices. Fig. 
EJb) shows the net result: a coarse-grained tensor net- 



FIG. 1. (a) Tensor network the ground state |T) of H on an 

infinite lattice. It is made of copies of tensor A and restricted 
to the upper half plane (:c,t + ), with a row of open indices 
at t = 0. (b) By coarse-graining the tensor network while 
leaving the open indices untouched, we obtain a new tensor 
network with tensors A' together with one row of disenta- 
glers and isometries, (c) Further coarse-graining of the tensor 
network produces new coarse-grained tensors A" and a sec¬ 
ond layer of disentanglers and isometries, (d) By iteration we 
obtain a full MERA approximation for state |T). 


work, with effective tensor A', together with a double 
row of disentanglers and isometries both for the incoming 
and outgoing indices. After 0(log 2 (/?)) iterations of the 
coarse-graining procedure, we obtain a MERA represen¬ 
tation of the thermal state, Fig. [2jc) made of 0(log 2 (/?)) 
double layers of disentanglers and isometries for both the 
incoming and outgoing indices, together with a central 
row of tensors. Disentanglers and isometries are isomet¬ 
ric tensors and thus do not affect the spectrum of eigen¬ 
values of the thermal MERA, which therefore depends 
exclusively on the central row of tensors [24]. 

The thermal MERA obtained from TNR resembles the 
form first suggested by Swingle in the context of holog¬ 
raphy [15], where a 1+1 conformal field theory is dual 
to a gravity theory in 3 space-time dimensions. In this 
context, the thermal MERA is interpreted as describing a 
space-like cross-section of a black hole space-time geome¬ 
try. A significant difference in our construction is the cen¬ 
tral row of tensors, which is absent in Swingle’s proposal 
0 and provides pp with the correct thermal spectrum 
of eigenvalues {e~^ Ei /Z}, where {Ei} are the eigenval¬ 
ues of H. This central row of tensors can be thought of 
as representing the Einstein-Rosen bridge connecting the 
two asymptotic AdS regions (25| , which seems to provide 
a manifestation of the ER=EPR conjecture [26|. 

Periodic chain of size L .— In our third construction, 
we use TNR to coarse-grain a tensor network for the 
ground state of Hamiltonian H on a periodic chain of size 
L, see Fig. [3](a) . The network consists of a semi-infinite, 
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(a) open boundary (b) 




FIG. 2. (a) Tensor network on an infinite strip of finite 

width /3, with two rows of open indices. It is proportional to 
the thermal state (b) By coarse-graining the tensor 

network while leaving the open indices untouched, we obtain 
a new tensor network with tensors A' together with an up¬ 
per and lower row of disentaglers and isometries, (c) Futher 
coarse-graining produces a thermal MERA. 


vertical cylinder of width L, with a row of open indices. 
After about 0(log 2 (L)) coarse-graining steps the size of 
the system has effectively become 0(1), see Fig. [3(c), and 
we have 0(log 2 (A)) layers of the MERA connected to a 
semi-infinite cylinder of 0(1) width. This semi-infinite 
cylinder can be understood as the infinite product of a 
transfer matrix T. The dominant eigenvector of T leads 
to the ground state of ff, whereas subdominant eigen¬ 
vectors describe low energy eigenstates. 

To illustrate the computational possibilities offered 
by the new algorithm, we consider the one-dimensional 
quantum Ising model with transverse magnetic field both 
at finite (3 for an infinite chain, and at zero temperature 
for a finite periodic chain of length L [27|. The calcu¬ 
lation required less than 5 minutes on a 2.5Ghz dual 
core laptop with 4Gb of memory (MERA bond dimension 
X = 10). First, for L = oc, Fig. [4](a) shows the expecta¬ 
tion value of the energy density ^thermal = tr(ppH)/L as 
a function of the inverse temperature /?, while Fig. 0(b) 
displays, at critical magnetic field, the crossover between 
polynomial decay of correlations at short distances (due 
to quantum fluctuations at criticality) and their exponen¬ 
tial decay at longer distances (due to finite temperature 
statistical fluctuations). Then, for /3 = oo, Figs. 0(c)- 
(d) show, respectively, the low energy spectra of ff as a 
function of the inverse system size 1/L, and the energy 
and momentum of low energy states for L = 1024. In all 
cases, an accurate approximation to the exact result is 
obtained. 

Discussion .— Since its proposal a decade ago Q, the 




FIG. 3. (a) Tensor network on a semi-infinite vertical cylin¬ 

der of finite width L and with a row of open indices, propor¬ 
tional to the ground state of if on a periodic chain made of L 
sites, (b) Result of coarse-graining the initial tensor network 
while not touching its open indices, (c) MERA connected to 
a semi-infinite vertical cylinder of 0(1) width. Inset: transfer 
matrix T of this cylinder. The eigenvectors of T with largest 
eigenvalues correspond to the low energy eigenstates of H. (d) 
MERA for the ground state/low energy excited states of if, 
where the top tensor is an eigenvector of the transfer matrix 
T. 


MERA has been regarded as a variational class of states. 
Accordingly, one attempts to approximate the ground 
state |T) of if by optimizing the variational parame¬ 
ters contained in the disentanglers and isometries of the 
MERA, for instance by iteratively minimizing the ex¬ 
pectation value of if [8, l9| . Such energy optimization is 
costly (it may require thousands of sweeps over scale) and 
prone to becoming trapped in local minima. Moreover, 
there is no guarantee that the end result is an approx¬ 
imation to the ground state |T) - one just has a wave- 
function with, hopefully, reasonably low energy. Here 
we have argued that, instead, an approximate MERA 
representation of the ground state |4/) can be obtained 
with TNR by coarse-graining an initial, quasi-exact ten¬ 
sor network representation of |4/) [29]. This only requires 
one sweep over scale, and it is therefore computationally 
much more efficient j30[. In addition, at each coarse- 
graining step TNR introduces a truncation error 0 that 
can be explicitly computed. If this truncation error is 
sufficiently small, then one can certify that the resulting 
MERA approximates the true ground state |4>) within 
that small error [31]. 


Importantly, TNR acts locally : the coarse-graining of 
the tensor network at point (x,r) only depends on the 
tensors in an immediate neighborhood 17]. Let us men¬ 
tion three consequences for the resulting MERA. (i) Since 
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FIG. 4. (a) Thermal energy per site (above the ground 

state energy) as a function of the inverse temperature (3, for 
the quantum Ising model H m JA XiX i+ 1 + A JA Zi in an 
infinite chain, for different values of magnetic field A. Con¬ 
tinuous lines correspond to the exact solution, (b) Connected 
two-point correlators at the critical magnetic field A = 1, as a 
function of the distance d , for several values of (3. Continuous 
lines correspond again to the exact solution, (c) Low energy 
eigenvalues of H for critical A = 1 as a function of 1/L. (d) 
Low energy spectrum of H for critical A = 1 and correspond¬ 
ing momentum (in unites of 2ir/L) for L = 1024 sites, which 
appear organized according to the conformal towers of the 
identity I (red), spin a (green), and energy density e (blue) 
primary fields of the Ising CFT, 28]. Discontinuous lines 
in (c) and (d) correspond to the finite-size CFT prediction, 
which ignores corrections of order L ~ 2 . 


TNR is not aware of the system size L or inverse temper¬ 
ature /?, it produces the same tensors A, A ', A ", • • • and 
disentanglers and isometries for the ground state |T) in 
the thermodynamic limit (L = f3 m oc) as it does for the 
states pp and |4/( L )) at finite /J or L. Thus, a single TNR 
calculation produces an accurate MERA approximation 
for all these states, (ii) In the absence of translation in¬ 
variance, where a different tensor Ai may be required for 
each site i of the lattice, the coarse-graining of different 
parts of the system can be conducted in parallel, leading 
to a massive reduction in computational time, (iii) At a 
conceptual level, locality of TNR implies the validity of 
the theory of minimal updates 33]. This theory, which 
asserts that only certain parts of the MERA need to be 
changed in order to account for a localized change in the 
Hamiltonian H 33], is particularly useful in the study of 
systems with boundaries, defects, or interfaces [9, 14J. 

In Summary, in this Letter we have shown that, when 
applied to the Euclidean path integral restricted to sev¬ 


eral geometries, TNR 0 produces a MERA for the 
ground state and thermal states of a quantum Hamil¬ 
tonian, and have explored a number of consequences of 
this result. We conclude by briefly mentioning two more 
implications. First, thanks to this connection, TNR in¬ 
herits from the MERA its ability to define an RG flow 
in the space of wave-functions and Hamiltonians |fj [§]. 
Notice that extracting the emergent physics from these 
RG flows should be easier (both conceptually and com¬ 
putationally) than extracting it from the RG flow in the 
more abstract space of tensors, Eq. [lj Second, although 
we have focused on quantum systems, TNR can also be 
applied to statistical partition functions 0- Therefore 
the present construction extends the MERA formalism 
(including strategies to extract universal critical proper¬ 
ties (lCHl4j) to classical statistical systems. 
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Appendix A.— A tensor network for the Euclidean 
path integral 

In this section we describe how to obtain a conve¬ 
nient, quasi-exact tensor network representation of the 
partition function/Euclidean time path integral Z = 
tr ( E~P h ) (or Euclidean time evolution operator e~P H ), 
starting from the short-ranged Hamiltonian H of a ID 
quantum system on the lattice. This tensor network is 
used in the main text as the starting point to obtain 
MERA representations for the ground state and thermal 
states of H in infinite and finite systems, see also Ap¬ 
pendix B. It is also used in Appendix C. 

For simplicity, we assume that H is a sum of nearest 
neighbor terms only, 

H = Y / Ki+i ( 2 ) 

i 

[Longer- (but finite-) range terms can be treated with a 
slightly more complicated scheme.] The first step is to 
split H into two contributions, H = H even + H 0 dd 

Heven = ^ ^ 1, L/odd = ^ ^ 1} (3) 

even i odd i 

and use a Suzuki-Trotter decomposition [34] to approx¬ 
imately express e@ H as the P-fold product of operators 

eleven anC [ e e #odd ? 

e^ H « (e“^even e -^oddj P ? 6 = ft/P < 1. (4) 

This introduces an error of order 0(/3e), which therefore 
vanishes in the limit of small e/large P. [One can ob¬ 
tain an error 0(/3(e) n ), n > 1, by using a higher order 
Suzuki-Trotter decomposition [33]]. Since H even is a sum 
of terms that act on different sites and therefore com¬ 
mute, e eHeven is simply a product of two-site gates, and 
similarly for e /3 ^ odd , 

g C-f^even __ g +1 g e -f^odd —— g +1 

even i odd i 

(5) 

Each two-site gate e i s a tensor made of four in¬ 

dices. Therefore we have obtained a tensor network rep¬ 
resentation of e~P H , see Fig. [5][a). Using singular value 
decompositions (in the vertical and horizontal direction, 
as indicated in the inset of Fig. [5j) we can then transform 
this tensor network into a new tensor network made of 
four-legged tensors A connected according to a square 
lattice pattern. 

However, by construction e~ ehi ^ i+1 is very close to the 
identity, e~ ehi ’ i+1 = I + O(e), which implies that the 
square tensor network is highly anisotropic. It is then 
convenient to coarse-grain in the Euclidean time direc¬ 
tion, until space and time direction have become qualita¬ 
tively equivalent, in the sense e.g. that the ordered sin¬ 
gular values of a tensor in both directions decay roughly 
in the same way. This is illustrated in Fig. [6j 



FIG. 5. (a) Tensor network for the Euclidean path integral 

obtain by multiplying small (Euclidean time) two-site gates of 
the form e ~ ehi ’ i+1 . (b) Same tensor network after decompos¬ 
ing each two-site gate using a singular value decomposition, 
according to the inset, (c) Final tensor network for the Eu¬ 
clidean path integral in terms of copies of a single tensor A. 



FIG. 6. (a) Tensor network for the Euclidean path integral 

produced in Fig. [5l made of copies of a single tensor A. (b) 
We can introduce isometries (obtained e.g. from the singular 
value decomposition of a vertical pair of tensors A) to define a 
new tensor A. (c) In the new tensor network, each row of ten¬ 
sors A accounts for two rows of tensors A, in such a way that 
the anisotropy between space and Euclidean time directions 
has been reduced. The initial square tensor network discussed 
in the main text is the result of applying this compression in 
the Euclidean time directions a few times. 
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Appendix B.— Four geometries for an Euclidean 
time evolution 

In this section we review the well-known fact that 
the Euclidean time evolution operator is pro¬ 

portional to the thermal state p ^, 



where L denotes the length of a periodic chain and /3 
is the inverse temperature. Here we use H^ and H = 
(j eno ^ e a Hamiltonian on a periodic lattice of 
size L and on an infinite lattice, respectively. 

We consider four geometries, obtained from combina¬ 
tions of infinite and finite L and /3. These geometries are 
represented in Fig. [71 The partition function Z'i L ^ corre¬ 
sponds to a torus of size L and p in the horizontal and 
vertical directions, respectively, which are parameterized 
by space and Euclidean time coordinates (x,r). Then 
e~PH (L) corresponds to the vertical cylinder obtained by 
cutting the torus at r = 0, as represented in Fig. [TJd). 

At infinite length L, this cylinder is replaced by an 
infinite strip of finite width /?, representing e ~^ H , which 
is proportional to the thermal state pp of Fig. 0b), with 



FIG. 7. Four geometries on which we define the tensor net¬ 
work for e ~ (3H . (a) Upper half plane (x,r + ), corresponding 
to the ground state |T) of if on an infinite lattice, (b) Hor¬ 
izontal strip of finite width /3, corresponding to the thermal 
state e ~^ H . (c) Semi-infinite vertical cylinder of finite width 
L, corresponding to the ground state |\|A jL )) of if on a peri¬ 
odic chain of length L. (d) Vertical cylinder of finite width L 
and finite height /3, corresponding to a thermal state p^ on 
a finite periodic chain. 


Pp = lim pf\ 

L—>• oo H 


(7) 


Appendix C.— Tensor Network Renormalization, 
step by step 


At infinite inverse temperature j3 (or zero temperature) 
we obtain a projector 

|*W><*W| = lim —-, (8) 

/3—>-oo y\^) 

onto the ground state of H^ L \ In other words, for 

any state |0o) such that (T^ L ^|0o) 7^ 0, 

|\I/( L )) oc lim e f3H( } |0o). (9) 

/3—>-oo 

Therefore, for finite L, the Euclidean time evolution op¬ 
erator on a semi-infinite vertical cylinder is proportional 
to the ground state of H^ L \ Fig. 0c). 

Finally, for infinite f3 and L, the semi-infinite cylinder 
above is replaced with the upper half plane and produces 
the ground state |T) of H , Fig. 0a). 

We have therefore attached a geometry to each of the 
ground states and thermal states |T), pp, |\I/( L )), and 
In the main text we consider explicitly the coarse- 
graining of the tensor network for the first three geome¬ 
tries in Fig. [71 that is (a) - (c), and the fourth one [the 
finite cylinder (d)] is left as an exercise that follows from 
the previous cases. 


In this section we review the detailed sequence of local 
replacements that define the coarse-graining transforma¬ 
tion of TNR. For convenience, we have introduced some 
cosmetic changes compared to Ref. [17 ]. 

The actual sequence is shown in Fig. [8j It 
involves a number of intermediate tensors (namely 
B,u,VL,VR,x,y,w) in addition to the final coarse¬ 
grained tensor A! . When applied to an infinite square 
tensor network made of copies of tensor A , Fig. [9j these 
local replacements produce a coarse-grained tensor net¬ 
work made of copies of a new tensor A', where there is a 
copy of A' for every four initial tensors A. None of the 
intermediate tensors appear in the coarse-grained tensor 
network. We can think of TNR as implementing a global 
change of scale, from (x,r) to (x',r') where 


-A x' = x/2, r —y t — t/2. 


( 10 ) 


Let us now consider instead applying TNR to coarse- 
grain the restriction of the initial square tensor network 
on the upper half plane, as represented in Fig. [TO] Recall 
that this restriction defines a projector onto the ground 
state of H. By applying the sequence as local replace¬ 
ments of Fig. 0 we see that most of the upper half plane 
ends up in a coarse-grained square tensor network again 
made of tensors A !. However, near the open boundary of 
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the tensor network first a row of disentanglers u , and then 
a row of isometries w, accumulate. These disentanglers 
and isometries cannot annihilate with Hermitian conju¬ 
gates vJ and because those would belong to the lower 
half plane, which is not available. 

As discussed in the main text and Fig.l of the main 
text, iteration of the coarse-graining produces a full 
MERA for the ground state |T) of H. Let us now give a 
geometrical interpretation of this construction. It turns 
out that TNR allows to implement local scale transfor¬ 
mations on the lattice [36], that is, a discrete version 
of certain conformal transformations in the continuum. 
From this perspective, the MERA results from (the dis¬ 
crete version) of the conformal trasformation that maps 
the upper half plane into the hyperbolic plane H 2 . 

If in the continuum we start with a flat metric in the 
upper half plane (x,r), the infinitesimal square line el¬ 
ement is given by dx 2 + dr 2 , which we can rewrite as 
r 2 (dx 2 /r 2 + dr 2 /r 2 ). Then a local change of scale (Weyl 
re-scaling) by r _1 (that is, by an amount that depends on 
the Euclidean time coordinate r) produces a new metric 

dx 2 + dr 2 (dx 2 /r 2 + dr 2 /r 2 ) (11) 

which corresponds to H 2 . Finally, we have also replaced 
the Euclidean time coordinate r with the scale coordi- 



FIG. 8. Sequence of elementary steps involved in TNR. (a) 
Each block of four tensors A is approximately replaced with 
an auxiliary tensor B, as well as disentanglers u and isometries 
vl and vr. The last three are determined through a local 
optimization [l7], that is, one that only knows about the four 
tensors A and is therefore insensitive to the rest of the tensor 
network, (b) Definition of the auxiliary tensor B in terms of 
A, u,vl, and vr. (c) Tensor B is approximately replaced by 
the product yx~ 1 . (d) Tensors x and y are such that y~ x Bx is 
equal to the identity, and are therefore also determined locally, 
(e) A combination of tensors is approximately replaced with 
tensor A' and isometries tc, ud. The isometry w is determined 
through a local optimization, (f) Definition of tensor A'. 



FIG. 9. (a) Square tensor network made of infinitely many 

copies of tensor A. (b) Plaquettes of four tensors A are re¬ 
placed by a tensor B and other tensors according to Fig. [8][a). 
(c) Pairs of disentanglers u and id annihilate according to 
md = I. (d) We decompose each tensor B according to Fig. 
He), (e) Tensor A' and isometries w are introduced accord¬ 
ing to Fig. H e )- (f) Pairs of isometries w and ltd annihilate 
according to icitd = /, producing the final tensor network, 
which contains only copies of the tensor A!. 


nate s = log 2 (r), which labels the different layers of the 
MERA. 

Returning to the discrete geometry of the square ten¬ 
sor network, notice that the net effect of iterating TNR 
on the (discrete) upper half plane as discussed above is 
that we are introducing a change of scale where locally 
the metric is changed by an amount proportional to the 
distance r = 2 s to the boundary, that is exponential in 
the scale variable s, as in the above conformal transfor¬ 
mation. 


Appendix D.— TNR versus energy minimization 


In this section we compare the accuracy and efficiency 
of two methods that allow us to obtain a MERA approxi¬ 
mation to the ground state |T) of a many-body Hamilto¬ 
nian H : (i) TNR and (ii) energy minimization mm- 
We use the ID quantum Ising model, with Hamiltonian 


h = XiX i+ i + a Zi, 


( 12 ) 


where X and Z are Pauli matrices, as a benchmark model 
for comparison. The ground state of the Ising Hamilto¬ 
nian H at critical magnetic field strength, A = 1, is rep¬ 
resented using a scale-invariant MERA obtained in these 
two ways. The scale invariant MERA considered here 
consists of three transitional layers followed by infinitely 
many copies of a scale-invariant layer |8|, 11. Il2j. 

TNR .— In the first approach, we start by generating 
the Euclidean path integral for iL, as described in sec- 
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FIG. 10. (a) Ground state |T) of if, represented by a tensor 

network on the upper half plane made of copies of tensor A. 
(b) Approximately equivalent tensor network obtained by re¬ 
placing blocks of four tensors A with the auxiliary tensor B 
and disentanglers and isometries, according to Fig. IH a), (c) 
The disentanglers annihilate pairwise according to uu^ = 7, 
everywhere except at the open boundary, (d) We approxi¬ 
mately replace tensor B by the product yx ~ 1 according to 
Fig. |8lfc). (e) An approximately equivalent tensor network 
is obtained by introducing tensor A' and the isometry w, 
according to Fig. [DJe). (f) The isometries annihilate pair¬ 
wise according to ww^ = /, everywhere except near the open 
boundary. The net result is thus a coarse-grained tensor net¬ 
work on the upper half plane with tensors A' and a row of 
disentanglers u and isometries w at the open boundary, which 
define a layer of the MERA. 


(a) (b) 



FIG. 11. (a) Relative error in the ground energy (per-spin) for 
the ground state of the ID quantum Ising model at criticality, 
comparing MERA of bond dimension x optimized either using 
TNR or variational energy minimization, (b) Relative error 
in the ground energy from (a), now plotted as a function 
of the leading order computational cost of the optimization 
method, which scales as 0(x 6 ) for TNR and 0(x 9 ) for energy 
minimization. 


tion A. Then we apply four coarse-graining iterations of 
TNR, after which the flow in the space of tensors has 
reached an approximate scale-invariant fixed-point Af p , 


see Eq. 1 in the main text. A scale-invariant MERA 
is then built from tensors (disentanglers and isometries) 
produced while implementing TNR, as described in the 
main text. Here the first three iterations of TNR generate 
the three transitional layers, while the fourth iteration of 
TNR is used to build the scale-invariant layer. 

Energy minimization .— In the second approach, we 
use the energy minimization algorithm of Ref. [8|, 11, uM 
to directly optimize a scale-invariant MERA (with three 
transitional layers preceding the scale-invariant layers) by 
iteratively minimizing the expectation value of the H. 

The comparison between the ground energy error SE 
obtained by the two methods is presented in Fig. [TlJ 
For an equivalent bond dimension x> energy minimiza¬ 
tion produces a MERA with smaller error SE than TNR, 
by a factor k « 10 that is roughly independent of x- That 
energy minimization would produce a more accurate es¬ 
timate of the ground state energy than TNR was to be 
expected, for two reasons. First, the energy minimization 
algorithm directly targets the ground energy as its figure 
of merits, whereas TNR does not. Second, the energy 
minimization algorithm optimizes the MERA by sweep¬ 
ing back and forth over length scales/layers many times, 
whereas TNR considers each length scale only once, with¬ 
out returning back to optimize shorter length scales after 
coarse-graining larger length scales. 

However the cost of optimization also differs between 
the approaches: the leading order cost of energy mini¬ 
mization scales as 0(x 9 ), whereas the cost of TNR scales 
as 0(x 6 )- There is also a significant difference in the 
number of optimization steps required to reach a con¬ 
verged result. The energy minimization algorithm re¬ 
quires many more steps to converge, as the algorithm 
sweeps over scales thousands of times, compared to TNR, 
which treats each scale only once. The difference between 
the two approaches in terms of the number of optimiza¬ 
tion steps required is difficult to quantify as it depends on 
many factors, but here we find that energy minimization 
requires roughly an order of magnitude more steps than 
TNR. Fig. [TTlf bl compares the energy error SE between 
energy minimization and TNR as a function of the lead¬ 
ing order computational costs of the approaches. Even 
though this comparison does not take into account the 
difference in the number of optimization steps between 
the approaches, TNR is already seen to be more compu¬ 
tationally efficient in producing a MERA approximation 
to the ground state within some specified accuracy SE 
than the energy minimization algorithm. 


Appendix E.— TRG yields a tree tensor network 
(TTN) 

The main result of this paper is to show that when 
TNR is applied to a tensor network representing the Eu¬ 
clidean path integral of a quantum Hamiltonian H (re- 
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stricted to the upper-half plane), it generates a MERA 
approximation for the ground state of H. In this sec¬ 
tion we show an analogous result for TRG: we show that 
when applied to the same tensor network, it generates 
instead a tree tensor network (TTN) approximation for 
the ground state of H. 

This result, which extends to several generalizations of 
TRG, offers an alternative perspective to understand the 
qualitative differences between TRG and TNR, and pro¬ 
vides a robust theoretical argument for why TRG can¬ 
not realize a critical fixed-point tensor while TNR can 
0. Indeed, the failure of TRG to produce a critical 
fixed-point tensor can be understood to follow from the 
well-known fact that the TTN ansatz (with fixed bond 
dimension) cannot reproduce the logarithmic scaling of 
entanglement entropy, Sl ~ log(L), found in the ground 
states of ID quantum systems at a critical point. 

For simplicity, we will derive this result not for the 
original TRG, but for an improved version known as 
higher order tensor renormalization group (HOTRG), 
which was proposed in Ref. [38] . However, the same con¬ 
clusions can also be shown to hold for the square-lattice 
form of the original TRG algorithm proposed in Ref. [5]. 
We begin with brief summary of HOTRG, see Fig. [T2j 

Consider a square lattice tensor network made of copies 
of a tensor A with components Aijki and bond dimension 
y. la HOTRG, one coarse-grains the tensor network first 
in one direction and then in the other. We shall begin 
with the horizontal axis. In the first step, left-right pairs 
of tensors A are multiplied together, obtaining the six- 
index tensor A , 

Ai,ji,j 2 ,k,h,i 2 = ( 13 ) 

p 

as shown in FigIT2fa). We now should truncate the dou¬ 
ble indices j = (jhjb) and l = (M 2 ) of A, which are of 
dimension y 2 , down to dimension y. In HOTRG this 
truncation is performed using a generalized form of the 
SVD called the higher order singular value decomposition 
(HOSVD), a method of decomposing an N rank tensor 
into the product of a core tensor and a set of isometric 
tensors. Using the HOSVD, the tensor A is approxi¬ 
mately decomposed into a core tensor B (with all four 
indices of dimension y), and a pair of isometries w and 
(with their three indices also of dimension y), 

x 

Ai,ji,j2,k,h,l2 ~ F W l’h,l2 B iJ,k,l( w jJi,h)*> ( 14 ) 

3 , 1 =1 

see also Fig. fl2l a). [Note that we are assuming, for 
simplicity, that the network is invariant under spatial re¬ 
flection along the horizontal axis such that the top and 
bottom isometries w are conjugates of one another, al¬ 
though the same result could also be reached without 
this simplifying assumption.] Isometries w satisfy the 


A A 
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FIG. 12. (a) A left-right pair of tensors A is multiplied 

to form a new tensor A, which is then decomposed using a 
truncated HOSVD into the product of a core tensor B and a 
pair of isometries w, see also Eqs. fl3l and IT4l (b) The prod¬ 
uct of an isometry w and its conjugate produces the identity, 
see Eq. 1161 (c) A up-down pair of tensors B is multiplied 
to form a new tensor D, which is then decomposed using a 
truncated HOSVD into the product of a core tensor A' and a 
pair of isometries v. (d) The product of an isometry v and its 
conjugate produces the identity, (e) The sequence of trans¬ 
formations from (a-d) above is used to coarse grain a 4 X 4 
network of tensors A, yielding a 2 x 2 network of tensors A' 
surrounded by isometries w and v. 


following relation, 

x 

w i,ji,j2( w jji,j2) ( 15 ) 

j 1,32*1 

as depicted in Figll2lbb After the HOSVD, we obtain 
a new square lattice tensor network composed of tensors 
B and rescaled by a factor of 2 in the horizontal direc¬ 
tion. Following the same steps but now in the vertical 
direction, up-down pairs of tensors B are multiplied to¬ 
gether to form a six-index tensor D, which can then be 
decomposed using the HOSVD into a product of a core 
tensor A' and a pair of isometries v and v 1, see Figure). 
[Again, for simplicity, we have assumed that the network 
is invariant under appropriate spatial reflection.] The 
combination of horizontal and vertical contraction steps 
has taken the initial square lattice network of tensors A 
to a new network of tensors A', one which is rescaled by 
a factor of 2 on both axes. This constitutes a single itera¬ 
tion of the HOTRG approach. By iteration, we generate 
a sequence of increasingly coarse-grained lattices defined 
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FIG. 13. (a) Square lattice tensor network representing 

the ground state of a quantum Hamiltonian if on a ring of 
L sites through an infinite Euclidean time evolution. The 
tensor network describes a semi-infinite vertical cylinder of 
width L and with a row of open indices, (b) The result of 
coarse-graining the initial tensor network with one iteration 
of HOTRG, see Fig. fl2f eb while not touching the L open 
indices at the bottom, (c) After coarse-graining with a sec¬ 
ond iteration of HOTRG, a TTN connected to a semi-infinite 
vertical cylinder of 0(1) width is produced. Inset: transfer 
matrix T of this cylinder. The eigenvectors of T with largest 
eigenvalues correspond to the low energy eigenstates of H. 
(d) TTN for the ground state/low energy excited states of if, 
where the top tensor is an eigenvector of the transfer matrix 
T. 


in terms of tensors, 




') ( w" ,v ") 

J A" ( 4 J . 


(16) 


Let us now apply HOTRG to a square lattice tensor 
network representation of the ground state of a Hamilto¬ 
nian if on a periodic chain of size L, as we did in the main 
text for TNR. The network consists of a semi-infinite, ver¬ 
tical cylinder of width L, with a row of open indices at the 
bottom, see Fig. fl3T a). Each coarse-graining iteration 


rescales the bulk of the network by a linear factor of 2, 
while leaving a row of isometries w along the open bound¬ 
ary, see Figll3l After about log 2 (L) coarse-graining steps 
the size of the system has effectively become 0(1), and 
we have 0(log 2 (L)) layers of the TTN connected to a 
semi-infinite cylinder of 0(1) width. This semi-infinite 
cylinder can be understood as the infinite product of a 
transfer matrix T ; taking the dominant eigenvector of T 
as the top tensor of the TTN leads to the ground state 
of H . Thus the HOTRG approach, applied to the Eu¬ 
clidean path integral of a quantum system iF, has been 
shown to generate a TTN approximation for its ground 
state. 

Since the introduction of TRG, a significant improve¬ 
ment in renormalization schemes for tensor networks has 
been the use of the so-called environment to achieve a 
more accurate truncation. The second renormalization 
group (SRG) proposed in Ref. [HI uses to knowledge from 
the environment, obtained approximately by iteratively 
sweeping over all length scales, to allow the tensor trun¬ 
cations to be performed so as to minimize the global er¬ 
ror, as opposed to the local error minimized in TRG and 
HOTRG. Ref. [HI introduced a version of SRG based 
upon HOSVD, called the HO SRG, which we will now 
consider. The use of the environment allows HO SRG 
to more accurately discern which degrees of freedom can 
be truncated, and which should be retained, at each RG 
step. This allows for a smarter choice of isometries. How¬ 
ever, HO SRG is structurally still based upon the same 
blocking scheme as HOTRG, Fig lEZT e). It follows that 
HOSRG also leads to a TTN representation of ground 
and excited states as depicted in Figll3j The only dif¬ 
ference being that the isometries w that compose the 
TTN have been better chosen with HOSRG than with 
HOTRG, as they minimize the global error. In particu¬ 
lar, for the same bond dimension HOSRG may pro¬ 
duce a more accurate TTN approximation to the ground 
state than HOTRG. 

However, as it happens with HOTRG (and TRG), 
HOSRG cannot generate a proper RG flow in the space 
of tensors - that is, one that at criticality flows to a criti¬ 
cal fixed-point tensor. Indeed, HOSRG produces a TTN, 
which is known to only be able to accurately describe a 
critical system if the bond dimension grows with increas¬ 
ing scale. But a bond dimension that grows with scale is 
incompatible with having a fixed-point tensor. 

























































